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Abstract. 

We present N-body/SPH simulations of the evolution of an isolated dwarf galaxy including a detailed model 
for the ISM, star formation and stellar feedback. Depending on the strength of the feedback, the modelled dwarf 
galaxy shows periodic or quasi-periodic bursts of star formation of moderate strength. The period of the variations 
is related to the dynamical timescale, of the order of 1.5 10* yr. We show that the results of these simulations are 
in good agreement with recent detailed observations of dwarf irregulars (dirr) and that the peculiar kinematic 
and morphological properties of these objects, as revealed by high resolution HI studies, are fully reproduced. We 
discuss these results in the context of recent surveys of dwarf galaxies and point out that if the star formation 
pattern of our model galaxy is typical for dwarf irregulars this could explain the scatter of observed properties 
of dwarf galaxies. Specifically, we show that the time sampled distribution of the ratio between the instanteneous 
star formation rate (SFR) and the mean SFR is similar to that distribution in observed sample of dwarf galaxies. 

Key words, methods: numerical, methods: n-body, galaxies: dwarf, galaxies: evolution, galaxies: ISM, galaxies: 
irregular 



1. Introduction 

The nature of the processes regulating star formation in ir- 
regular galaxies is poorly understood. Whereas there is at 
least some understanding of star formation in regular spi- 
ral galaxies, this is less so for irregulars. For spiral galax- 
ies the guiding observation that the star formation rate 
(SFR) is related to the gas surface density by the Schmidt 
law has given rise to a number of competing theories that 
reproduce the general feat ures o f star formation i n larg e 
spiral galaxies (Elmegreen I2OO2', Dopita & Rvder ll994 . 
These systems seem to be regulated by large scale gravi- 
tational instabilities. Star formation in irregular galaxies 
has proven to be more difficult to understand. Irregulars 
have a widely varying SFR, spanning 4 ord ers of magni- 
tude for the normalized SFR/area fHunter [l997|) . possi- 
bly due to the fact that gas thermodynamics, governed 
by varying heating a nd co oling; processes, plays the deci- 
sive role fElmegreen |2002() . But why do some irregulars 
have very high SFRs relative to their mass, while others 
hardly show any activity? Are there any intrinsic proper- 
ties of the galaxies that can explain this disparity between 
SFRs or do all dwarf galaxies exhibit episodes of high star 
formation? 

In recent years a number of studies have highlighted 
these questions by investigating samples of dwarf galaxies 
and comparing their properties as derived from photom- 
etry, HI and Ha observations. Van Zee (j200CL ■2001.) in- 



vestigated a sample of isolated dwarf galaxies and found 
no strong correlation between star formati on an d inde- 
pendent physical parameters. Hunter et al. l|l998j) tested 
different regulating processes, amongst which disk insta- 
bilities, thermal and shear regulated star formation, but 
found that none could explain patterns of star formation. 
Stil (1999) investigated the relation between star forma- 
tion and HI gas kinematics. 

On the other hand, detailed studies of a number of 
nearby dwarf galaxies have fully revealed the complex 
structure of the ISM in these systems. H igh re solution 
aperture synthesis mapping (e.g. Kim et al. 1998, Wilcots 



& Miller 1998, Puche et al. 1992. Walter & Brinks 20011) 
of their HI has shown the interstellar medium (ISM) of 
these dwarfs to be a frothy structure, with holes of varying 
sizes, shells and filaments, even extending far beyond the 
optic al rad ius. From velocity dispersion studies (Young 
et al. Il997^ the presence of cold and warm neutral com- 
ponen ts pre dicted by the two phase model for the ISM 
rPield 119651) has been deduced. Comparison with Ha and 
UV observations shows that the dense walls of these holes 
are the sites of star formation (Walter et al. l200l|) . and 
sugg est 'ch ains' of successive star forming sites (Stewart 
et al. l2000|) . The cause of the holes in the HI distribution 
seems to be the energy input from ionizing radiatio n, stel- 
lar winds and supernovae, although Rhode et al. lll999^ 
and Efremov et al. l|l998|l discuss other possible scenar- 



ios. 
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Together these two types of observations have painted 
a picture of the complex interaction between star for- 
mation and the ISM of these systems that is formidable 
to capture theoretically. Some early attempts have been 
made to understand star formation qualitatively by the 
application of stochastic self propagating star formation 
(SSPSF ) mod els to dwarf galaxies (Gerola & Seiden .lQSOL 
Comins flQsi . While they probably capture some general 
characteristics of star formation, they are phenomenologi- 
cal and do not include the underlying physics of the ISM. 

Efforts to investigate the influence of star formation on 
the ISM of dwarf galaxies have mainly concentrated on the 
effects of large (central) bursts and on questions concern- 
ing the ejection of gas an d the distr ibutio n of metals(e.g. 
Mac-Low & Ferrara Tl999l Mori et al. ll997j) . Recently there 
also have been some simulations adressi ng the question of 
survival of small galaxies (Mori et al. |2002() . Generally 
these simulations have not tried to set up a self consistent 
model for the ISM and star formation, but prescribed a 
certain SFR. 

The importance of a good model for the ISM and 
feedback has been recognized by a number of authors. 
Andersen & Burkcrt ( 2000|) formulated an extensive model 
for the ISM in terms of a phenomenological model for 
the interstellar clouds. Their model showed self regula- 
tion of the SFR and the y foun d moderate fluctuations in 
SFR. Berczik & Hensler (2003^ incorporated such a cloud 
model into a chemodynamical galaxy evolution code. 
Semelin and Combes ( 20o3) formulated a model with sim- 
ilar characteristics, representing clouds by sticky particles, 
but did not apply these to dwarf galaxies. Springel & 
Hernquist ((2003) formulated a subgrid model for the mul- 
tiphase interstellar medium, producing a quiescent self- 
regulating ISM. However, relatively little effort has been 
directed towards resolving the normal evolution of irreg- 
ular dwarfs and providing the connection with detailed 
studies of single systems and extensive unbiased samples. 
Nevertheless dwarf galaxies are good test systems for ex- 
ploring star formation in galaxies: they are dynamically 
simple systems in the sense that they do not exhibit spi- 
ral density waves or shear. Furthermore their small size 
means that simulations can follow the various physical 
processes at finer linear and density scales. As the small 
scale physics of star formation and feedback presumably 
do not differ between normal and dwarf galaxies, we can 
use results obtained from these simulations and apply the 
same methods to larger systems. Here we present results of 
a simulation of the evolution of a normal dwarf irregular 
galaxy including a detailed model for the ISM, star for- 
mation and feedback. The distinguishing characteristics 
of this work are that the model for the ISM we employ 
does not explicitly postulate the presence of a two phase 
medium, rather it forms it as a result of the physics of the 
model. Furthermore we take special care in formulating 
a star formation model that is solely based on the Jeans 
instability, and we formulate a feedback scheme that gives 
us unambiguous control over the strength of the feedback. 
We will discuss the results of the simulation both in rela- 



tion to detailed observations of comparable single systems, 
as well as in the context of recent surveys of dwarf galax- 
ies. 

2. Method 

We employ an N-body/SPH code for the evolution of a 
general astrophysical fluid on ga lactic scales, extended 
from TreeSPH (Hernquist & Katzlis?), to simulate the 
evolution of an isolated dwarf galaxy. We use the conser- 
vative SPH formulation of Springel & Hernquist l|2002|) . 
Main features of our code are: a realistic model for the 
ISM solving for the ionization and thermal balance for the 
neutral and ionized components of the ISM, star forma- 
tion based on a gravitational instability model for clouds, 
and a new method of including feedback for SPH. We will 
summarize the features of the code with an emphasis on 
the aspects most relevant for the present work. 

2.1. Model for the ISM 

Our model for the ISM is, although simplified, qualita- 
tively similar to the model for the Cold Neutral Medium 
(CNM) and Warm Neutral medium (WNM) of Wolfire et 
al. (,1995. .20031) . We consider a gas with arbitrary but fixed 
chemical abundances Xi, scaled to the target met allicit y 
from the solar abundances of Grevesse & Sauval l)l998j) . 
We solve for the ionization and thermal evolution of the 
gas. The various processes included are given in Tabled A 
similar model to that employ ed her e was used by Gerritsen 
& Icke l(l997j) and Bottema ()2003|) for galaxy simulations. 
The main differences are the following: we use more ac- 
curate cooling, that is calculated in accordance with the 
chemical composition, we have included a solver for the 
ionization balance, and we use the full p hotoel ectric heat- 
ing effic iency as given in Wolfire et al. l|l99!4l . Gerritsen 
& Icke l|l997l) found that the structure of the resulting 
ISM depended strongly on the ionization fraction they as- 
sumed, as this strongly influences the cooling. We do not 
have to assume an ionization fraction, as we calculate it 
(on the other hand, we do assume a cosmic ray ionization 
rate that is poorly constrained). Our use of the full heat- 
ing efficiency means that FUV heating will become less 
efficient for high radiation fields, due to grain charging. 
In our model supernova (SN) heating is more important 
in re gulati ng star formation than it was for Gerritsen & 
Icke lll 997|) . 

A concise overview of the ISM model is given in Fig. ^ 
The plots in this figure show that as density varies, the 
equilibrium state of the gas changes from a high temper- 
ature/high ionization state (T = 10k K, Xe ~ 0.1) at 
low densities, to a low temperature/low ionization state 
(T < 100 K, Xc < 10^^) at high densities. In between is a 
density domain where the negative slope of the P-n rela- 
tion indicates that the gas is unstable to isoba ric pre ssure 
variations, the classic thermal instability fField ll965|) .The 
shape of these curves and hence the exact densities of the 
thermal instability vary locally throughout the simulation 
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H, He only; primary 
& secondary ionizations 


1 


Collisional 


H, He only 


3 



Radiative recombination H, He only 3 

CIE assumed for metals 

Table 1. Overview of the processes included in the ISM model used. For H and He ionization equihbrium is explicitly 
calculated, for other elements collisional ionization equilibrium (CIE) is assumed. Both the heating and cooling strongly 
depei id on the ionization fractionate. Exact expressioi is ado pted for the vario us pro cesses can be found in: 1) Wolfire 
et al. ll995i 2) Raga et al. ll997l . 3) Verner & Ferland ll99d 4) Silva & Viegas 1200X1 



according to the conditions of UV and supernova heat- 
ing. The gas in the simulation may be out of equilibrium, 
although the timescales for reaching equilibrium are gen- 
erally short, < 10^ yr. In principle the cooling proper- 
ties of the gas depend on the local chemical composition. 
In practice, only small metall icity g radients are observed 
in dwarfs (Pagel & Edmunds 1 1 98 I L so we take constant 
metallicity (Z = 0.2 Z©). Potentially more serious is the 
fact that we assume a constant cosmic ray ionization rate 
C throughout the galaxy. The low energy (« 100 MeV) 
cosmic rays that are important for heating and ioniza- 
tion of the ISM, have relatively short mean free paths( 
< 10 pc) so probably C will vary substantially across a 
galaxy. However, the exact sources, let alone production 
rates, of these low energy cosmic rays are not well known. 
Hence a satisfying model for the distribution of cosmic 
rays is difficult to formulate. We take a 'standard' value 



of C = 1.8 10 



-17 



For this value cosmic ray heating 



only becomes important in the outer parts of the galax y 
(it may be substantially too low, see McCall et al. '2003^. 

The FUV luminosities of the stellar particles, which 
are needed to calculate the local FUV field used in the 
photo electric heating, are derived from Bruzual & Chariot 
l|1993l and updated) population synthesis models for a 
Salpeter IMF with cutoffs at 0.1 Mgand 100 M©. In the 
present work we do not account for dust extinction of UV 
light. 



2.2. Star formation and Feedback 

We use t he star formation recipe of Gerritsen & 
Icke lll997l) . A region is considered unstable to star forma- 
tion if the local Jeans mass Mj is smaller than the mass 
of a typical molecular cloud Mj-of ~ IO^Mq. The rate of 
star formation is set to scale with the local free fall time: 
f,ftff = 

10 seems reasonable 



, The delay factor fgf is uncertain, 

but from observations a value fgf 



(Zuckerman & Palmer ri974|) . Once a gas particle is de- 
termined to be forming stars, a fraction Sgf — 0.25 of the 
mass is converted to stars. This sets a minimum to the 
star formation efficiency (a neighbourhood of gas parti- 
cles of about 64 particles that has become capable of star 
formation will thus form at least a fraction .25/64 — .004 
of stars). The actual efficiency of star formation is de- 
termined by the number of stars needed to quench star 
formation locally by the UV and SN heating and is deter- 
mined by the cooling properties of the gas and the energy 
input from the stars. 

The recipe is certainly not unique, we have based our 
choice on the following considerations: it is simple and 
based on presence of substructure and the driving role of 
gas self gravity, and it reproduces the Sch midt l aw with- 
out actually imposing it( Gerritsen & Icke 1199*^ . An ad- 
ditional advantage is the fact that the Jeans mass crite- 
rion for star formation also prevents the simulation from 
violating the resolutio n requ irements for self g ravitating 
SPH (Bate & Burkert Il997i Whit wort h [l99a). It does 
assume substructure to be present (actually its assump- 
tion is more restrictive still, namely that substructure is 
mainly at Mref sized clouds, but this is not overly impor- 
tant). For this type of simulation one is always restricted 
by the limited ability to follow the starforming process, so 
we are forced to adopt a phenomenological description at 
some level. Another concern that could be adressed is that 
the SF according to this recipe is independent of metallic- 
ity, other than that induced by the metallicity dependence 
of the cooling. Furthermore, we assume the stellar IMF to 
be universal. 

The inclusion of feedback in SPH simulations iscom- 
plicated by the so called overcooling problemfKatz [l99^ . 
A number of methods has been devised to overcome this, 
based on the retur n of energ y as e ither thermal (Thacker 
& Couchman "2000", Gerritsenll997() or kinetic( Navarro & 
White .1993. Springel & Hernauist l2003i) . None has been 
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Fig. 1. Overview of the ISM model: Equilibrium plots of 
(a) temperature T, (b) electron fraction and (c) pres- 
sure P as a function of density n for two different values 
of the UV field Go, drawn line: Go = 1, dotted: Go = 10 
( xl.6 10~^ ergs cm~^ s~^ ) 

entirely satisfactory as the former requires the artificial 
suppression of cooling, while the latter suffers from too 
strong effects for feedback, unless an unrealistically low SN 
efficiency is adopted or the affected particles are decoupled 
from the gasdynamics. Here, we employ a new method. It 
is based on the creation at the location of young stel- 
lar clusters of a zero mass SPH particle. Its contribution 
to gas forces on neighbouring gas particles is determined 
by taking the formal zero mass limit with constant parti- 



cle energy Esn of the usual SPH equations of motion(see 
Pelupessy et al. l2004l for more details). The amount of 
feedback( the energy EgN given to the zero mass SPH 
particle) is an important and uncertain parameter; we es- 
timate it by multiplying the number of type II supernovae 
per mass of stars formed, O.OO9M0~"'^, with the effective 
energy per supernova, « lO^^erg, thus assuming that 90% 
of the energy is radiated away, a value which comes from 
more detailed simulations of the effect of su pernova and 
stellar winds on the ISM (Sihch et al. llQQfil) . and is also 
used in ot her sim ulations of galaxy evol ution f e.g. Semelin 
& Comb es Eool Springel & Hernauist Eool Buonomo et 
This energy is then gradually added to the su- 
pernova particle over a period of 30 Myr, the lifetime of SN 
progenitors, after which it is ceases to exist. This model 
for the mechanical energy input from stellar clust ers re- 
sembles the one adopted by Oey & Clarke l|l997|) . who 
have shown that it is fully consistent with the differen- 
tial size distribution of superbubbles in a.o. the SMC and 
Holmberg II. We do not include the effect of SN la, which 
are mainly important for the chemical enrichment. 

2.3. Initial conditions 

Here we will restrict ourselves to follow the evolution of 
an isolated dwarf galaxy in an attempt to understand the 
interplay between star formation and the ISM. These pro- 
cesses play a role for dwarf galaxies in general, and indeed 
in all places were star formation is happening. It is impor- 
tant to keep in mind, however, that dwarf galaxy evolu- 
tion is easily influenced by external factors such as infall 
of gas, collisions, and recently receiving much attention, 
ram pressure stripping by the inter galactic medium and 
tidal stripp ing by galactic potent ials ( see e.g. Marcolini 
et al. l2003l Mori fc Burkert Eoool and Mayer et al. l200li 
Pasetto et al. l2003l for recent work on these). Nevertheless 
the isolated case remains very relevant as it allows a 
straightforward comparison with suitably selected sam- 
ples of o bserve d galaxies which exist in the literature(e.g. 
van Zee l200ll) , unaffected by the uncertainties of addi- 
tional external influences. The Understanding of the com- 
plex processes of star formation and feedback gained from 
the isolated case can then be aplied to more general evo- 
lution scenarios. 

We thus setup a model dwarf galaxy resembling cur- 
rent dirrs. Although these exhibit a wide range of mor- 
phologies, they are very similar in their averaged radial 
profiles to scaled down versions of ordinary disk galaxies. 
Hence we take for the initial condition a three component 
model for a small disc galaxy, consisting of a gas disk, 
a stellar disk and a dark halo, modelled loosely on the 
properties of DDO 47 and similar dIrrs. The gas disk we 
construct with a radial surface density profile, 

S = I]g/(l + R/Rg), (1) 

with central density Sg — 0.01 Mo/kpc^ and radial scale 
Rg — .75 kpc, truncated at 6 kpc. It is somewhat involved, 
due to the thickness of the gas disk, to solve its vertical 
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structure exactly. Hence, before we start the simulation 
proper, we set up the galaxy with an quadratically rising 
gas scaleheight and we it run isothermally(T = 10* K) and 
without star formation and feedback for some time until 
it settles in equihbrium. This results in a scaleheight set 
by hydrostatic equilibrium of hg = 100 pc at the center to 
hg = 1 kpc at R = 6 kpc.An exponential stellar disk, 



Pdisk(R., z) 



2h, 



■ exp(— R/Rd)sech^(z/hz] 



(2) 



with central surface density cto — 0.3 M0/kpc^, Rd = 
0.5 kpc and vertical scale height h^ = 0.2 kpc, is con- 
structed as in Kuijken & Dubinski ((iQQfTl) . The total mass 
of the gas is Mg = 2 x IO^Mq and the total mass of the 
stellar disk is Md = 1.5 x IO^Mq. Both the gas disk and 
stellar disk are represented by N = 10^ particles. The ages 
of the initial population of stars are distributed according 
to a constant SFR of 0.007 Mg/yr. 

Dwarf galaxies are amongst the most dark matter dom- 
inated objects, exhibiting dark to luminous matter ratios 
of 10-100. Their rotation curves are best fit by fl at cen - 
tral density cores ( Flores & Primack ll994L Burkert ll995|) . 
Therefore we take a halo profile 



Phaio(r) = Po 



exp(-r2/r^) 



1 



(3) 



with core radius 7 = 2 kpc, cutoff radius rc = 20 kpc 
and central density po = 2 x 10^ Mo/kpc"^, for a total 
mass of Mhaio = 15 x lO^M© and a peak rotation velocity 
of about 50 km /s. The profile Q is very similar to the 
Burkert l|l995fl profile for dwarf galaxies, differing in its 
asymptotic behavior for r — > oo. This will only give signif- 
icant deviations well outside the region of interest for our 
simulations. We represent the dark halo by a static po- 
tential. This is deemed sufficient for the dynamical mod- 
elling presented here, as we evolve the galaxy in isolation 
and the perturbations in the gaseous and stellar disk are 
expected to have only minor impact on the halo struc- 
ture. Furthermore, discreteness noise of a particle halo 
can induce bars to form, and excite spiral structure or 
bending modes in the stellar and gaseous disk( Kuijken & 
Dubinski . 199,5. Hernauist .1993.) . 

3. Simulation Results 

In Fig. 12 a range of snapshots of the simulated dwarf 
galaxy is shown, taken after about 1 Gyr of evolution 
and spaced about 50 million years apart. Plotted are UBV 
composite images with HI contours overplotted. If we look 
at the contours and maps of HI of the upper panels of 
Fig. Owe see that the HI distribution of the galaxy is dom- 
inated by holes of varying sizes. Dense shells of HI combine 
to form big HI clouds, which correspond to the sites of 
intense star formation. Sometimes structures resembling 
spiral arms form, but generally there are no spiral density 
waves in gas or stars, and, due to the low shear, a floc- 
culent spiral structure also does not develop. As the cold 



clouds move about in the galaxy they are in turn destroyed 
by mechanical feedback from the stars to perpetuate the 
cycle of star formation, cloud destruction and formation. 
Clearly the ISM of the galaxy is in a very dynamical state. 
This can only be fully appreciated while looking at the 
time sequence of the complete simulation. The ISM of 
dwarf irregulars is continously stirred and material from 
different radii is mixed through the action of supernovae 
and stellar winds. While this may have been expected be- 
forehand from energy considerations, it is nonetheless an 
important and often overlooked fact. As such it should 
not be surprising that irregulars often show very little ra- 
dial metallicity gradients (apart from the fact that metals 
may be lost from hot outfiows). Note that the presence of 
the hydrodynamics of the ISM is a fundamental difference 
between classical SSPSF and phcnomological models and 
our full dynamical simulations. 

If we look at the UBV images we see that the appear- 
ance of the galaxy varies with time. Episodes of strong 
star formation( e.g. at 1024 Myr) are followed by quies- 
cent phases( 1071 and 1118 Myr ). During star formation 
bursts the galaxy is dominated by a few very active star 
formation sites, during quiet times the galaxy fades and 
seems relatively featureless. If we compare the stellar and 
HI distribution of Fig.|21to observed systems, it is striking 
how similar these are. For example Holmberg II and NGC 
42 14, wh ich have been studied in great detail (Puche et 
al. ll992l Walter & Brinks .2001.) . show an HI distribution 
with a similar number of holes of similar sizes. Also the 
extent of structure in the gas outside the optical disk is 
reproduced in the simulations. Large regions of star forma- 
tion are associated with the highest contours of HI (corre- 
sponding to a column density of about N = 2 x lO^^cm^^), 
which is also found in comparisons of optical and HI mor- 
phology for observed systems. 

As we have a complete representation of the neutral 
phases of the ISM we can also compare our simulation 
with the HI kinematics of observed systems. In Fig. Owe 
have plotted maps of the HI density and the (line of sight) 
velocity dispersion of the gas. The gas typically has ran- 
dom velocities of the order of 3-10 km/s, with regions of 
higher (up to 20 km/s) velocity dispersion associated with 
expanding bubbles, in good agreement with observations 
(^Stil ll999i Young et al. ll997l) . 

The evolution of the SFR is shown in Fig. 0] It shows 
a gray scale plot of the temporal evolution of the (az- 
imuthally averaged) star formation density and below 
it the resulting total SFR. We see that after an initial 
transient period the galaxy settles in a mode of period- 
ically varying star formation, with a minimum SFR of 
about O.OO3M0/yr while the peak SFR is about a fac- 
tor of 10 higher. Also plotted are the results for a run 
with reduced feedback strength. In this case variations are 
smaller in amplitude, and quasi-periodic. Even lower feed- 
back strengths would probably quench the variations com- 
pletely, but in that case morphology and kinematics of the 
gas would no longer match observed systems. In the upper 
panel of Fig. ^we see the pattern of propagating star for- 



6 



Pelupessy, van der Werf, Icke; Bursts of SF in Irregulars 





Pelupessy, van der Werf, Icke: Bursts of SF in Irregulars 
2.0 




(kpc) 



(kpc) 



Fig. 3. HI surface density & velocity dispersion map in quiescent phase(left, taken at 1165 Myr) and in burst phase( 
right, at 1211 Myr). 



mation. The azimuthaUy averaged density plot may give 
the false impression that star formation is happening in in- 
ward and outward moving rings; however visual inspection 
shows that it moves around in patches and partial rings: if 
some region starts to form stars, nearby dense regions will 
be triggered to form stars, often moving the star formation 
in a particular direction along a dense filament or bridge. 
Furthermore we see that periodic bursts happen only rel- 
atively close to the centre of the galaxy, whereas in the 
outer parts stars form at a more or less constant lo w rate . 
This behaviour is expected (see Ehlerova & Palousf^^) 
because star formation can only be triggered in the ex- 
panding shells around the holes induced by feedback if 
the column density is high enough, N « io^°~^^ cm~^. So 
we see that the galaxy has three different SF regimes: an 
outer region where the gas density is generally not high 
enough for the CNM to form (beyond the so called thermal 
cut off, see Elmegreen & Parravano 1994 and Gerritsen & 
Icke Il997ll so SF will be strongly supressed there; a re- 



gion at intermediate radii where the gas density is high 
enough to be thermally unstable, and where star forma- 
tion proceeds at a steady pace; and a central region that is 
both thermally unstable and unstable to shell instabilities, 
where the dominant mode of SF is triggered star forma- 
tion. Fig.Elshows a plot of the resulting radial dependence 
of the star forming density, clearly visible are the three dif- 
ferent regions with different star forming behaviour. The 
bursting region is confined within 1.5 kpc, while star for- 
mation extends further out, stopping abrubtly at « 4kpc, 
well short of the edge of the gasdisk. Interestingly, De Blok 
& Walter ( 2003) found evidence for similar low level star 
formation outside the optical disk of NGC 6822. 

The cyclical pattern of star formation raises a number 
of further questions: What is the driver of these variations? 
What determines the period and amplitude of the oscilla- 
tion? Why is star formation apparently synchronized over 
the whole galaxy? While some previous models of galax- 
ies have shown variations of the star formation rate, these 
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Fig. 4. Star Formation History of a simulated dwarf galaxy. Upper panel: density plot of the azimutally averaged star 
formation rate. Lower Panel: total star formation rate. The dotted line indicates the SFR for a rmi with 50% reduced 
feedback strength. 
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Fig. 5. Radial profile of the mean star formation density( 
drawn line, dotted lines are exponentials with scalelengths 
of .25 and 1. kpc). 



were either connected to an imposed time scale (the du- 
ration of supernova feedback, Wada & Norman l200lh or 
du e to ra ndom scatter of a probabilistic model (Gerola et 
Significant quasi r andom variations in SF rate 
were found by Carraro et al. (■2001) to happen during for- 
mation of dwarf ellipticals (dEs), in their case caused by 
local depletion of the starforming gas reservoir, happen- 
ing on timescales of « 10^~^ yr. Similar variations driven 
mai nly b y internal proccesses were found by Pasetto et 
al. ll200.1^ even for dirrs subjected to tidal stirring. In our 
case the timescale of the variations is related to a phys- 
ical time scale: the period of the oscillation is about 165 



Myr, which is slightly larger than double the dynamical 
time scale of the darkhalo, (tdyn — \J itop ^ Myr). 
This suggests that supernova feedback provides the kick 
for the oscillation, expelling gas and increasing the veloc- 
ity dispersion and thus inflating the gas disk, after which 
the gas falls back for the next cycle on a dynamical time 
scale. The period is thus determined by the density of the 
halo (the main mass component of dwarf galaxies), while 
the amplitude is related to the strength of the supernova 
feedback. Note that this would not be sufficient to ac- 
count for the regularity of the variations: in general one 
would expect that gas from different parts of the gas disk 
would experience kicks of different strength, and would 
fall back on different time scales, giving rise to a more 
chaotic pattern of star formation. However dwarf galaxies 
are in solid body rotation and thus their potential cor- 
responds to a harmonic oscillator, which means that ma- 
terial flung out of the disc falls back in the same time, 
irrespective of location and speed, synchronizing the SFR 
over the whole galaxy. This effect certainly adds to the 
phenomenon shown in Fig. 01 namely that the overall star 
formation cycle looks remarkably regular. 

4. Star formation: comparison with observations 

When we compare the star formation in ou r model and 
observed SFR (for example in Van Zee l200ll) we see that 
the agreement is quite good: our values fall well within 
the range of observed rates of O.OOlMo/yr to O.lMo/i/r, 
although admittedly the star formation rate in our model 
does depend fairly strongly on poorly constrained quanti- 
ties such as the shape of the IMF, the amount of energy 
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released in SN and stellar winds, and the various param- 
eters regulating star formation. It is however encouraging 
that we can use reasonable values for the parameters and 
get SFR that are close to observed rates. 

In our model the star formation rate is seen to vary 
over time. This is difficult to verify in real dwarf galaxies. 
One would have to derive accurate SFH over hundreds of 
Myr for galaxies, and variations of within a factor sa 3 of 
the mean as observed in our model would not be too con- 
spicous in the data. For most dwarf galaxies the long term 
star formation histories determined from global colours 
indicate an approximately constant SFR. Detailed exam- 
inations of nearby resolved dwarfs do show evidence o f 
variations in the SFR (Dohm-Palmer et al. Il998l |2002() 
of the correct order of magnitude. Strong additional evi- 
dence for variations in the SFR is the observed scatter in 
the ratio of current star formation to past star formation: 
the Scalo b parameter. If the variations in our model were 
generic for isolated dwarf galaxies, this would explain the 
observed distribution of b. To see this we have plotted in 
Fig. El the histo grams of SFR/(SFR)t for our simulation 
together with the distribution of b for a sample of iso- 
lated dwarf galaxies as derived by Van Zee ( 200l|. Clearly 
these two are qualitatively the same. There are two main 
differences: (1) the distribution of our model extends to 
higher and lower ratios and (2) the galaxy sample contains 
galaxies with 6 > 4, which are extreme bursting systems. 
The first diff'erence would be remedied by decreasing the 
strength of feedback. We also have plotted the results for 
a run with 50% less feedback. In this case the amplitude 
of variations is somewhat too small(although both are in 
good agreement with the observed distribution given the 
uncertainties inherent in deriving it, see discussion in Van 
Zee i2001il . For the second point we should mention that 
our simulation really only samples a limited number of 
cycles; the bursts may be rarer events that have not (yet) 
happened in the simulation. Ofcourse it also possible that 
high powered bursts are due to outside influences, such as 
interactions or infalling gas clouds, ingredients of galaxy 
evolution not present in our model. 

The distribution of is a highly degenerate indicator 
for the SF history, and we are aware that the above com- 
parison presupposes a number of properties for the popu- 
lation of dwarf galaxies: the assumption that the dwarf 
galaxies form a homogeneous population; variations in 
SFR are assumed to be of similar strength and shape, 
independent of galaxy properties ( although the b distribu- 
tion is independent of the period); and that the variations 
in SFR are periodic. Nevertheless we think that this is 
the simplest explanation for the observed scatter and that 
our model captures a number of essential features of the 
variations in star formation given the characteristic shape 
of the observed b distribution, namely: moderate burst- 
ing, long periods of low SF, and the 'peakiness' of the 
SF enhancements. It is also difficult to see what would 
induce non-periodic variations in isolated galaxies, espe- 
cially since no properties independent from the SFR that 
correlate with the SFR (Van Zee^QQli) have been found. 



0,3 



0/ 




12 3 4 ^5 

Fig. 6. distribution of SFR/(SFR)t for a sample of iso- 
lated dwarfs (squares with error bars), a simulation with 
full feedback (drawn line) and a simulation with 50% f eed- 
back( dotted line). Data is taken from Van Zee l200ll 



5. Discussion &. Conclusions 

The model we employed to calculate the evolution of a 
single dwarf galaxy is quite extensive and complicated 
and necessarily contains a number of free parameters, 
some of which we tried to constrain using observational 
data or theory ( as in the case of the supernova energy), 
yet some remain only constrained by physical intuition. 
Therefore, it is worthwhile to summarize again the fea- 
tures of observed dwarf systems successfully reproduced 
by the model: 1) the morphology and kinematics of the 
HI distribution, 2) the two phase structure of the ISM, 2) 
the spatial pattern of star formation, 3) the star forma- 
tion rate, 4) consistency with the observed distribution of 
Scalo b parameter. These provide strong and independent 
checks of our modelling approach as they result from the 
intricate interplay of the model for the ISM we employed, 
the scheme for star formation and the method of returning 
mechanical energy to the ISM. These successes give some 
confidence that the model captures the essentials of dwarf 
galaxy evolution and to its predictive power. 

Comparing our model with recent simulations of dwarf 
galaxy evolution ( e.g. Berzci k et al. l2003[ Ma yer et 
al. l20nil Andersen & Burkert 1200(1 P asetto et al. '2003', 
Mori et al. [1997, 1999, Carraro et al. l200l|) we see that 
our work differs from previous work in two crucial aspects: 
1) the modelling of the physics of the neutral ISM, 2) the 
treatment of supernova feedback. We follow the evolution 
of the WNM and CNM of the ISM exphcitly, this process 
being the upper part of a cascade leading down to star 
formation. Star formation in our model happens in cold( 
T <; 200 K) and dense (n = 1 - 10 cm^^) gas. These are 
realistic sites for star formation as we know that molecu- 
lar clouds are embedded in neutral envelopes. The further 
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stages of star formation, molecular cloud formation and 
collapse, are not included as they require prohibitive reso- 
lution and the inclusion of additional physics. These pro- 
cesses are only implicit in our star formation recipe, but 
at least our methods follows collapsing gas to the point 
that it has experienced a transition to a cold, dense state 
from which it may be trusted to form stars with rates and 
efficiencies that are constrained by observations. Some au- 
thors have tried to bypass this problem by formulating a 
phenomenological ISM model in terms of 'sticky' particles 
representing molecular clouds (Andersen & Burkcrt 2000), 
sometimes in addition to a smooth SPH c ompo nent rep- 
resenting the WNM (Semelin & Combes [20021 Berczik 
et al. |2003). They also succesfully reproduce a selfreg- 
ulated ISM, including some effects, like evaporation of 
molecular clouds that are probably not well represented 
in our model. In these models stars form from the cloud 
particles, imposing a Schmidt law, not, as in our model, 
from the consideration of the the instabilities in the ISM. 
Furthermore, our model has a consistent representation of 
the ISM linking phases by physical processes, rather than 
prescriptions. 

Simulations of dwarf galaxy evolution that include su- 
pernova feedback have typically included this as a local 
heating term. This method for implementing supernova 
feedback is not effective in forming the structures associ- 
ated with stellar feedback and simulations using it do not 
show any effect of SN feedback. This is a well understood 
numerical artefact and some authors have devised meth- 
ods to prevent the radiative loss of mechanical energ y that 
is the root of the problem (Springel & Herna uist 1200,'^ 
Thacker & Couchman kOOd Gerritsen & Icke Il997t the 
problem is a lso no t exclusive to SPH type codes see e.g. 
Fragile et al. l2003|) . They have not applied their methods 
to the evolution of dirrs so a direct comparison is not yet 
possible. However, the notion that the energy input from 
stellar winds and SN of young stellar clusters is crucial 
for the understanding the structure and kinematics of the 
ISM of dIrrs is borne out by our simulations: if feedback is 
not included, the ISM stays in a smooth disk with very low 
random motions( 2-3 km/s). The self propagating mode 
of star formation is absent in that case. The dwarf galaxy 
will only poorly resemble a real dirr. 

Chemical enrichment is not yet included in our code. 
We can estimate the importance of a changing chem- 
ical composition for our simulation. The total amount 
of metals produced during the simulation is Mz ~ 
0.015AM* = 1.6 X IO^Mq. This will raise the metallicity 
of the galaxy, under the assumption that the metals will 
be well mixed, with at most 25%, which is not entirely 
insignificant, although to first approximation the thermal 
evolution is independent of metallicity as both the cooling 
and UV heating scale with Z. Note also that metals may 
be lost from the galaxy in galactic winds. It does mean 
that we should include chemical enrichment to follow the 
evolution for longer timescales and especially if we want 
to investigate extremely metal poor galaxies. In principle 
the inclusion of chemical enrichment will also add further 



constrains comparing with observations, very important 
to asses the long term evolution of the model. Both these 
points, however, do not alter the conclusion for simula- 
tions as presented here. 

The methods we use can be applied more general. 
Questions that we plan to address are for example the 
following: 

— In view of the recent realization that there is a deficit of 
observed small galaxies as compared to predicti ons o f 
ACDM models of galaxy formation (Klypin et al. ll999() 
it is interesting to consider what happens for galaxies 
of progressively smaller mass. We are in the process 
of running a grid of models exploring this question, 
but some effects may be clear from the preceding: for 
smaller galaxies variations will be on longer time scales 
(scaling as P/7aL) ^^'^ bigger in amplitude (because the 
feedback becomes relatively stronger as the escape ve- 
locity decreases) . Ultimately halos will be too small to 
retain their ISM, leaving them bare. The mass ranges 
for which this happens, the influence of other galac- 
tic parameters and the timescales on which these pro- 
cesses take place will be of interest to validate cosmo- 
logical and galaxy formation models. 

— We also plan to investigate a possible relation between 
dwarf irregulars and dw arf ell ipticals. Alth ough r ecent 
simulations (Mayer et al. l200ll Pasetto et al. 2003) have 
shown for satellite galaxies that a transition from dIrr 
to dE or dwarf Spheroidal (dSph) is plausible as a re- 
sult of the action of tidal fields, it has not been conclu- 
sively determined whether a transition from dIrr to dE 
is to be expected in general. Various people ha ve put 
forward arguments in favo ur (Da vies & Phillipps'l988') 
and against (Bothun et al. 1198(1 Marlowe et al. .199&) 
such a descendancy for the dE. Although the results 
presented here do not immediately illuminate this 
question, we think that additional simulations of the 
sort presented here, testing a wider range of galactic 
properties and following the evolution for longer time 
scales may answer whether this is a viable scenario 
and whether we can put the various classes of small 
galaxies into an unified evolutionary framework. 

In summary, our model suggests that it is possible that 
a large part of the current dIrr population is in a quasi- 
periodic burst mode of star formation. The scatter in ob- 
served properties in this picture is mainly due to the fact 
that galaxies are observed at different phases of their evo- 
lution. The main difference between our model and classi- 
cal SSPSF models is that in our models variations are due 
to the interplay of stellar feedback and gasdynamics, the 
galaxy being periodically stirred by bursts of star forma- 
tion after which a quiescent period occurs during which 
gas falls back to the disk. In classical SSPSF models the 
variations are due to the fact that star formation is de- 
scribed by a correlated stochastic process and the small 
size of the galaxies, which induces large statistical varia- 
tions. Our model predicts that the amplitude of the vari- 
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ations depends on the strength of the feedback and that 
the period depends on the dynamical time scale. 
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